Heat conduction in one-dimensional Yukawa chains 
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Heat conduction in one-dimensional Yukawa chains is investigated. It is shown numerically that 
it has the abnormal heat conduction which is proportional to the system size. Effects of asymmetric 
external potential, the modified Frenkel-Kontorova one, on the heat conduction of system are also 
studied. It is found that the asymmetric property of external potential can induce asymmetric 
thermal conductivity that can be used to be a effective thermal rectifier. In certain of system 
parameters the heat flux are significantly different for two opposite direction. One can parametrically 
control the heat flux through this system by changing the potential strength and width. 
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Recent years many interesting works have been ad- 
dressed to the problem of heat conduction PJ 0, IE El 
IE IE IE 13 It is well known that heat conductivity 
of one dimensional (ID) lattice is an important problem 
which is related to microscopic foundation of the Fourier 
law. It also provides an abundant example for the mi- 
croscopic origin of the macroscopic irreversibility phe- 
nomenon. Detailed review of the problem is presented in 
the recent paper 

Except of extensive numerical studies focusing on the 
validity of Fourier law J = —ndT/dx for many mod- 
els such as the Fermi-Pasta-Ulam (FPU) model 0, , 
the ding-a-ling model H H|, the Frenkel-Kontorova (FK) 
model |E an d the Lorenz-gas model 0] and so on, 
some other different and important problems have also 
attracted people attention. For example most recently 
Casati et al. 0] investigated the possibility to control the 
energy transport inside a nonlinear ID chain connecting 
with two thermostats at different temperatures. They 
emphasized that controlling heat conduction by nonlin- 
earity opens a new possibility to design a thermal recti- 
fier, i.e., a lattice that carries heat preferentially in one 
direction. 

In this Letter we investigate heat conduction problem 
for a new type of model wherein there is a Yukawa- 
Coulomb interaction potential. It arises from dust- 
grain plasma which is extensively studied recently due 
to its new features with dust crystal or collective waves 
EE IH 13 III El EH Our study shows that the heat 
conductivity is abnormal in absence of external field as 
k oc N as in harmonic chains. Under FK field, however, 
it is normal which is similar to the cases of former studies 
but with harmonic or/and FPU interactions Surpris- 
ingly for the dust chains with real wake potential [l2l |. 
modified FK field, the heat conduction exhibits a mix- 



ture property of normal and abnormal behaviors. It is 
found that the asymmetric property of external wake po- 
tential can induce asymmetric thermal conductivity that 
can be used to be a thermal rectifier. In a certain of sys- 
tem parameters, by changing the strength and width of 
wake potential, the heat fluxes for two opposite direction 
can be different remarkably. 

In a general form the Hamiltonian of the considered 
system is written as 
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+ V{x l _ 1 x l ) + U(x i ) (1) 



where V(xi—iXi) stands for the interaction potential of 
the nearest-neighbor particles and U(xi) is a external on- 
site potential. For dust particles system we have the 
interaction potential (refer to EH) 
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where a is the parameter that ranges from to 1, Ad is 
the plasma Debye length and q is the dust grain charges, 
while r.ij is the distance between particles z-th and j-th. 
For simplicity we have assumed that the masses m<i and 
charges of all the grains q are the same and constant, 
while for convenience we shall choose parameter a that 
makes the neighbor dust particle is located at the bot- 
tom of V (see more clear in the following) . Note that the 
plasma electrons and ions do not appear explicitly in the 
simulations but their effects, which result in an attrac- 
tive force to dust, are taken into account in the effective 
interaction potential as the second term of Eq.(|2J). On 
the other hand in dusty plasma we know that there exist 
wake potential acting the dusts along downstream E2 
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where q t is the upstream test dust particle charge, M > 1 
is the Mach number of ion flow and ru is the distance of 
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FIG. 1: Temperature profiles for different external potentials 
of field-free (a), FK (b) and FK-modified (c). The imposed 
temperature are T+ = 0.12 and T_ = 0.08. The parameters 
are chosen as 71 = 1 and 72 = 0.3 (b), 2 (c). Averages are 
carried over a time interval 10 7 by dropping a transient 10 6 . 
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FIG. 2: Scaling of the heat flux JN with the number of parti- 
cles N in the Yukawa interaction chains for different external 
potentials of field-free (a), FK (b) and FK-modified (c). Sys- 
tem parameters and others are same as in Fig.l 



i-th particle from the test particle. The potential © is 
named as the FK-modified one in this paper because it 
is different from FK potential with varying amplitude. 

Let us turn to the side of the motion of system. A 
general approach is employed here that are used by many 
authors like Lepri et al. and Hu et al., namely, two Nose- 
Hoover thermostats are put on the first and last particle, 
keeping the temperature at T + and T_ , respectively. For 
simplicity and convenience, firstly we need to normalize 
the basic physical quantities: the mass by the particle 
mass rrid, the space distance by t he dust inter-particle 
distance ad and the time by the \J ma^/q 2 respectively. 
Then the equations of motion for particles are 

. 2 

xx = -£+xi + fi - / 2 , C+ = Tfr - 1 



'ii= fi-fi+h i = 2,...,N-l (4) 

. 2 

xn = -£,-xn + f N - Jn-i, £- = -pjr- - 1 

where /, = — V (x, — — U (a;,) is the force. 

Usually the fixed boundary conditions xq = const and 
xjv + i = const are assumed. Now the interaction poten- 
tial in normalized form can be written 

Vixi-^Xi) = (exp[-x(a; l - aji-i)] - a) /(af» - a^-i) (5) 

where parameter \ = a c i/Xd is the ratio of inter-particle 
distance to the plasma Debye length. Many studies show 
that when x is large the nearest neighbor interaction is 
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FIG. 3: Dependence of Jl/Jr, the heat flux ratio of opposite 
direction for the Yukawa- Coulomb model, on parameter 72, 
the strength of the FK-modified on-site external potential. 
The imposed temperature are T+ = 0.12 and T_ = 0.08 for 
Jr and T+ = 0.08 and T_ = 0.12 for J L . 71 = 1 is given. 

appropriate instead of particle-pair interaction |l5j . In 
our following numerical studies x = 5 is generally used. 
The particles are arranged initially at equilibrium po- 
sitions of V and the normalized inter-distance equals 
to 1, therefore, it means that parameter a, adjusting 
the attraction of plasma to dust particle, must satisfy 
a = (1 + x) ex P(~ x)- For the external FK-modified po- 
tential Eq.|j3J) we have its normalized form 



U(xi) =72 cos(7iXi)/xj, 



(6) 



where 71 = x/VM 2 - 1 and 72 = 2/3/(1 - M~ 2 ) with 
P — Qt/q the ratio of upstream test particle charge 
to downstream particles charge. Now we denote x^ as 
the position of i-th dust particle in the chain from the 
origin of test dust particle xot = 0. In correspon- 
dence for the normalized external FK potential one has 
U(xi) = 72 cos(7iXi). The different temperature profiles 
for different external potentials of field-free, FK and FK- 
modified are shown in Fig.l. The imposed temperature 
are T+ = 0.12 and T_ = 0.08 and 71 = 1. The potential 
strengths are 72 = 0.3 and 2 for FK and FK-modified, 
respectively. Obviously the temperature gradients are 
constructed globally only in the FK case. 

Usually the heat flux at the i-th position is given by 

Ji(t) = -n{Xi + Xi + i)fi + \. Numerically the time average 
J = (Ji(t)) is independent of the index for long time 
enough. We plot in Fig. 2 the heat flux JN vs the sys- 
tem size N for three different cases of field- free (a) , FK 
(b) and FK-modified (c). Obviously the abnormal heat 
conductivity k cx N is obtained for purely Yukawa inter- 
action chains, which is similar to the case of harmonic 
situation. In the Yukawa chains for standard FK exter- 
nal potential, although it seems unclear about whether 
there exists phase transition of heat conduction between 
normal and abnormal |g, |7j, our results indicates that 
the normal heat conductivity exists, which support the 
conclusion of normal heat conduction in Ref. |fjj. From 
Fig. 2 (c), however, we can see that it is abnormal glob- 
ally on heat conduction in the Yukawa chains under FK- 
modified on-site external potential while the local tem- 
perature gradient appeared in the left region of system, 
see Fig. 1(c). The dependence of Jl/Jr, the heat flux ra- 
tio of opposite direction for the Yukawa-Coulomb model, 
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FIG. 4: Momentum excitations for different types of external 

potentials of field- free [(a), (b)], FK [(c), (d)] and FK-modified 
[(e), (f)]. Pl(0) and Pr(0) denote the initial momentum kick 
on the left 1-th particle and right 100-th particle, respectively. 
For (c)-(f) the parameters of 71 = 1 and 72 = 2 are given. 



on parameter 72, the strength of the FK-modified on-site 
external potential is plotted in Fig. 3, where the power 
law Jl/Jr oc io - 24 '*' 2 is fitted in our numerical parame- 
ter regime. Our numerical experiments shows the valid- 
ity of this power law until to 72 = 4 which indicates that 
the difference of two directional heat current can possibly 
over one order. However in practice of dusty plasmas it 
can only be realized by a relatively small 72 due to the 
test dust charge equals or less than the charge of par- 
ticles in chains. For example, for given \ — 5, 71 = 1 
and /3 = q%/q ~ 1 it corresponds to the practice ion flow 
Mach number M = \/26 (usually in the sheath M > 1), 
and 72 ~ 2 which indicates Jl/Jr ~ 3. 

In order to see how the single kicked particle trans- 
fer energy to the others, in Fig. 4 we plotted the mo- 
mentum excitations with different external potentials of 
three cases of field-free, FK and FK-modified. We can 
see that in the case of field-free there exists a solitary- 
like wave which can carry energy independently and 
not transform it to others. In case of FK there exist 
phonons which can transfer heat to the others. And in 
the case of FK-modified there seems a combination of two 
kind waves mentioned above: solitary-like and scatter 
phonons. Moreover it seems that there are always more 
solitary-like waves which propagates in the left-direction 
than that in the right direction. For example, there are 



4 




FIG. 5: Dependence of the heat flux ratio of two cases, for the 
Yukawa chains with different inverse width of the FK-modified 
on-site external potential, on temperature T. 72 = 2 is given. 



three solitons in Fig. 3(f) but only one in Fig. 3(e). Maybe 
this is one reason why in this case the left-directional heat 
flux, Jl, is larger than that of right directional, Jr. 

On the other hand according to the expression of heat 
flux, the ratio of them of two-direction is approximately 




T + r > 1 \ 
expx(r> -r<) 

T- r < 



(7) 



where r> and r < are the average inter-particle distance in 
the region where the temperature is almost constant for 
[T+,T_], right-direction heat-flow J R , and [T_,T+], left- 
direction heat-flow, respectively. In general it is shown 
that r> is a little larger than r< by numerical experi- 
ments. It is not surprising that the more phonons are 
scattered in the case of [T + ,T_] than that of vise versa. 
Therefore even if r> is just a very little greater than r<, 
Jl/Jr would be greatly larger than one. For example, 
r> w 1.08 and r< ss 0.92, then we have Jl/Jr ~ 3. How- 
ever in the case of Harmonic interaction potential because 
of its symmetric quadratic property about equilibrium 
position the ratio of opposite directional heat flux is in 
order of ^/T + /T_, for example by choosing T + = 0.12, 



T_ = 0.08, the ratio of Jl/Jr ~ 1.23. Our numerical re- 
sults confirm this simple approximation analysis. Indeed 
if we instead of harmonic interaction chains the heat flux 
difference for two direction can not be over 23% under 
asymmetric potential. 

Finally we calculate heat flux in different potential 
width 7-f 1 (see Fig. 5). It shows that we can also con- 
trol the energy flow by changing the potential width 
7f 1 = \JM 2 — l/x through adjusting the ion flow veloc- 
ity, which is associated with M or and the plasma par- 
ticles density which is associated with \- Obviously the 
great difference of heat flux at low temperature caused 
by the width of potential smear out as the temperature 
increases to over 0.5. It is noted that the ratio of heat 
flux Jr(7i = 1)/Jr(7i = 7t) is about over two orders 
difference. 

In summary we have studied the heat conduction prob- 
lem in the system of ID Yukawa interaction chains with 
asymmetric on-site modified FK wake potential. We can 
parametrically control the heat flux through this system 
by changing the potential strength and width. Moreover 
it is also possible to adjust the heat flux by changing the 
potential from FK to modified FK in the Yukawa chains, 
thus we provide a simple method to change the prop- 
erties of the system, from a normal conductor obeying 
Fourier law, down to an almost perfect insulator. Some 
theoretical consideration has also been attempted to un- 
derstand for our results and it is found that in the case 
of Yukawa interaction potential, the asymmetric prop- 
erty about equilibrium position has played a crucial role 
to improve the ratio of opposite directional heat flux. It 
seems that there exist some modes which are combina- 
tions of solitary-like and scattering phonons. Certainly 
more formidable theoretical works is still challenging our 
research on this problem in the future. 
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